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Abstract Gauge theories with compact symmetry groups possess topologically non-trivial con¬ 
figurations of gauge field. This has dramatic implications for the vacuum structure of Quantum 
Chromo-Dynamics (QCD) and for the behavior of QCD plasma, as well as for condensed matter 
systems with chiral quasiparticles. I review the current status of the problem with an emphasis 
on the interplay of chirality with a background magnetic field, and on the observable manifes¬ 
tations of topology in heavy ion collisions, Dirac semimetals, neutron stars, and in the Early 
Universe. 
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1 Introduction 

Quantum Chromodynamics (QCD) possesses a compact gauge group SU(3) that 
allows for topologically nontrivial solutions for the gluon field. The existence of 
these solutions changes the vacuum structure - a superposition of an infinite set 
of topologically distinct states connected by tunneling instanton transitions ([Tj) 
becomes the “0-vacuum” of the theory (2,3). It is likely that topological effects 
are linked to confinement and chiral symmetry breaking in QCD, see (4) for a 
review. 
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Extended topological configurations of gluon and quark fields in QCD can be 
effectively probed by a background Abelian magnetic field. The interplay of QCD 
topology with the topology of zero modes of chiral fermions in a magnetic field 
leads to a number of surprising novel phenomena that are reviewed below; see the 
volume © and Refs. (6[[7|[8 ! 9) for complementary reviews of topics not covered 
here in detail. 

Under extreme conditions of high temperature and/or high baryon density, the 
vacuum of QCD changes its properties, and deconfinement and chiral symmetry 
restoration take place. This deconfinement transition is accompanied by the 
rapid change in the rate and nature of topological transitions connecting different 
topological sectors, as discussed in section [2| The heavy ion program reviewed 


recently in (10) opens a possibility to study these phenomena in experiment. 
Moreover, since the colliding ions create strong magnetic fields 0(10 m^) (11) 


(see (12) for review), the interplay of QCD topology with an Abelian magnetic 
background can be studied experimentally, as discussed in section [5j The exis¬ 
tence of strong magnetic fields in neutron stars allows to study novel topological 


phenomena in cold and dense QCD matter that are discussed in section 7.1 The 
opposite limit of high temperature and small baryon density is reached in the 


Early Universe, see section 7.2 


The recent discovery of Dirac semimetals - see (13) for a review — opens a 
fascinating possibility to study the topological phenomena involving 3D chiral 
fermions in condensed matter experiments. In section [6] we describe the recent 
observation of the chiral magnetic effect in ZrTes (14). 
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2 Above the barriers: Chern-Simons number diffusion 


2.1 Topology across the QCD phase transition 

The existence of the degenerate vacuum states with different Chern-Simons num¬ 
bers induces a P- and CP-odd “0-term” in the lagrangean of SU(N) gauge theory: 


where 


£qcd = -\G^{x)G a ^{x) - 9 q(x), 


q ( x ) = ^2 9 2 G* {x)G ailu {x) 


(1) 


( 2 ) 


m 


is the density of topological Chern-Pontryagin charge, and 6 angle is a new pa¬ 
rameter of the theory <©• The absence of observed global P and CP violation i 
QCD indicates that 9 is very close to zero - however the requirement of periodic¬ 
ity with respect to this angular variable puts important constraints on the QCD 
dynamics. 

At zero temperature, reconciling large N (number of colors) scaling with peri¬ 
odicity dictates (15) that the vacuum energy E{9) is a multi-valued function of 
9/N : 

E(9) = N 2 rnirife E > ( 3 ) 

where k is an integer. As 9 varies adiabatically, the ground state moves from one 
branch of the potential Eq. ^ to the other; note that this implies the existence 
of unstable vacuum states. 

At small 9 , we can write the free energy at a finite temperature T as 


F(6,T) - F(0,T) = - X (T) 9 2 f(6,T), 


( 4 ) 


where x(T) is the topological susceptibility. Expanding around 9 = 0, we can 
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further write the function f(9. T ) as 


f(6,T) = i + b 2 e 2 + ..- 


( 5 ) 


only even powers of 9 are allowed by the P- and CP-invariance of the ground state. 
The temperature dependences of topological susceptibility and of the moment b 2 
encode important information about the changes in the topological structure 
of the vacuum, and can be studied in lattice QCD simulations. The behavior of 
topological susceptibility as a function of temperature has been computed recently 
in SU(N) gauge theory; the results are shown in Fig. [I] One can see 
that the topological susceptibility rapidly decreases around the deconfinement 
transition, and this dependence sharpens at large N in accord with the scenario 


proposed in (18). The rapid decrease of topological susceptibility implies the 


existence of metastable P- and CP-odd vacuum states near T c (18). 


on 


To gain an additional insight, let us examine the lattice results 
the temperature dependence of the moment b 2 presented in Fig. [2j At high 
temperatures, color screening suppresses instantons with size /)> 1/T Ji9|20|2lt , 
so the dilute instanton gas picture should be appropriate, see for review. In 
this picture, the 9 dependence of the free energy is given by 


F(d,T) - F(0,T) ~ X (T) (l-cos0), (6) 

which yields a definite prediction for the moment b 2 . As can be seen from Fig. 
[2j the dilute instanton gas approximation is indeed appropriate for high tem¬ 
peratures in the deconhned phase; this regime is approached faster at large N. 
However at low T and around the deconfinement transition the topological struc¬ 
ture of QCD matter is not captured by the instanton gas picture - this implies 
that the effects of confinement on topological fluctuations are substantial and 
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have to be understood. 


2.2 Sphaleron rate in hot QCD plasma: weak coupling 


At finite temperature, the transitions between the vacuum states with differ¬ 
ent topological Chern-Simons numbers can be induced by a classical thermal 


activation process, so-called ’’sphaleron” (22). In electroweak theory sphaleron 


transitions cause the baryon number violation and may be responsible for at 


least a part of the observed baryon asymmetry in the Universe (23); for a review, 


see (24). In QCD plasma, sphalerons are abundant (25) and induce the quark 


chirality non-conservation. 

Unlike for the tunneling instanton transitions, the rate of the sphaleron tran¬ 
sitions r is not exponentially suppressed at weak coupling g, and in Yang-Mills 


theory with N colors is proportional to (26 27,28) 


T = const x ( g 2 N) 5 \n(l/g 2 N ) T 4 , 


( 7 ) 


with a numerically large pre-factor (29,30). Sphalerons describe a random walk 


in the topological Chern-Simons number Ncs'- 

(N CS (T) - N cs m 2 


T = lim 

T—>oc 


VT 


( 8 ) 


so that in a volume V and after a (sufficiently long) time period T we get the 
topological number ( Nq S ) = T V T. 

The expression Eq.Q can be qualitatively understood in the following way: 
since the sphaleron at the peak of the barrier separating Chern-Simons sectors 
is a purely magnetic field configuration, the factor ( g 2 T ) 3 is the inverse mag¬ 
netic screening length that determines the characteristic inverse volume of the 
sphaleron, and g 4 T log(l/ g) is the typical inverse time scale of the process. 
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The existing weak coupling computations are based on semiclassical field theory 
approaches, and do not provide a reliable result for the values of strong coupling 
a s ~ 1/2 relevant for quark gluon plasma produced in heavy ion collisions. A 


desperate extrapolation” (30) of the weak coupling results yields an estimate 


T ~ 30 a 4 s T 4 . 


( 9 ) 


2.3 Chern-Simons diffusion in non-Abelian plasma at strong cou¬ 
pling: holography 

At strong coupling, a valuable insight into the dynamics of non-Abelian confor¬ 


mal theories is offered by the holographic AdS/CFT correspondence (31,32,33). 


The holographic correspondence establishes duality between a strongly coupled 
conformal field theory in the boundary Minkowski space-time and classical super¬ 
gravity in the AdSs bulk space. The emergence of a black hole in the AdSs space 
corresponds to the formation of plasma with temperature equal to the Hawking 
temperature of the black hole. 

This approach allows to evaluate the Chern-Simons diffusion rate in conformal 
J\f = 4 maximally super-symmetric Yang-Mills theory (A f = 4 SYM) by consid¬ 
ering the bulk propagation of axions (excited by the topological charge operator 


on the boundary) in the background of the black hole (34). The absorption of 
the axions by the black hole with a cross section equal to the area of the black 


hole delivers the following result (34): 


T = 


(ff 2 AQ 2 

256 7r 3 


( 10 ) 


showing that topological transitions become more frequent at strong coupling, 
even though the dependence on the coupling is weaker than suggested by Eq0 











This means that Chern-Simons diffusion is not necessarily associated with semi- 
classical field configurations that are well defined only at small coupling and are 
usually washed out by quantum effects at strong coupling. Note that the large 
N behavior is the same in the weak and strong coupling limits ( N°). 

Of course, QCD is not a conformal theory, and the quark-gluon plasma at 
temperatures T < 2 T c significantly deviates from the scale-invariant behavior, 
as signaled by the expectation value of the trace of the energy-momentum tensor 


(0^)t = e — 3 P measured on the lattice (35) - the energy density e and pressure 
P of relativistic conformal plasma are related by e = 3 P. It is thus interesting 
to establish the effect of this deviation from conformality on the Chern-Simons 
diffusion rate. 


This question has been addressed (36) recently within the “improved holo¬ 


graphic QCD” (37) - a bottom-up approach describing the running QCD cou¬ 
pling and confinement by a properly chosen dilaton potential in the bulk. The 


resulting Chern-Simons diffusion rate (36) is significantly larger than that given 


by Eq.(10): at a s = 0.5, the authors find T(T C ) > 1.64 Tf at the transition tem¬ 


perature. This is about 36 times larger rate than given by Eq.(10) at a s = 0.5, 
and is about the same (within a factor of two) as the weak coupling extrapolation 

Eq.Q. 

It is important to identify configurations responsible for Chern-Simons diffusion 
at strong coupling. An explicit example of a topological solution in semiclassical 
AdS§ x S 5 supergravity (dual to the large N, strongly coupled J\f = 4 SYM) is pro¬ 
vided by the D-instanton (38,39) that is obtained from solitonic D-brane solutions 
by wrapping them around an appropriate compact manifold. D-instanton can be 


viewed as an Einstein-Rosen wormhole connecting two asymptotically Euclidean 
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regions of space-time, with the Ramond-Ramond (R<g>R) charge flowing down the 


throat of the wormhole (38). It describes a process of violation of the conser¬ 


vation of a global charge in the boundary theory (38). D-instantons have been 


proposed (40) as a source of multiparticle production in high energy collisions 


in strongly coupled A/"=4 SYM, extending on the earlier ideas of (41,42,43,44) 
based on the weak coupling approach. 


2.4 QCD phase diagram in the (T, /i 5 ) plane 

Because the chiral anomaly relates the net chirality of fermions to the topology 
of gauge fields, Chern-Simons diffusion generates the quark chirality imbalance. 
This imbalance can be quantified by the chiral chemical potential ^ = (/r# — 
Hl }/2 that describes an asymmetry between the Fermi energies of left- and right- 
handed quarks. The phase diagram of QCD in the (T, ^ 5 ) plane can thus shed 
light on the connection between confinement and topology. 

Moreover, since the chiral chemical potential does not lead to the fermion 


determinant sign problem (45), this phase diagram can be studied in the first 


principle simulations on the lattice. Note however that because of the chiral 
anomaly, the chiral chemical potential does not correspond to a conserved charge 
- the decay of chirality into gauge fields with non-trivial topology induces the 


so-called “chiral magnetic instability”, see discussion in section 7.1 


It has been argued in (|46,[47j 48]) that at finite fj ,5 the critical temperature 
decreases, and at some value of ^ the phase transition becomes first order via 


the second order end-point, in analogy to the behavior expected (49) in the 


(T, //) plane, see Fig. [ 3 ] The first lattice studies at finite ^ have recently been 


performed in SU{2) theory (50) with quark masses corresponding to the pion mass 
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in the range m v ~ 300-^500 MeV. Contrary to the model expectations (46, 47p8 ), 
the lattice results indicate that the critical temperature for the deconfinement and 
chiral symmetry restoration increases at finite pL §. It remains to be seen if this 
conclusion persists at physical quark masses - if it does, the result would imply 
a mechanism of interplay between confinement and topology not captured by 
current models. 


3 QCD in an external magnetic field 


3.1 QCD phase diagram in magnetic field 

An external magnetic field provides a unique way to probe QCD topology; it is 
thus instructive to study the phase diagram of QCD in a magnetic background. 
One of the key phenomena affecting chiral fermions in magnetic held is the “mag¬ 
netic catalysis” (51), see 0 for a recent review. This term refers to the enhance¬ 
ment of dynamical symmetry breaking induced by an external magnetic held. In 
a strong magnetic held, this phenomenon can be understood in terms of the ef¬ 
fective D —> D — 2 dimensional reduction. In particular, magnetic held promotes 
the formation of scalar chiral condensate - since the spins of fermions and an¬ 
tifermions on the lowest Landau level point in opposite directions, the formation 
of a spin-zero condensate is favored. Note however that the enhancement of sym¬ 
metry breaking in magnetic held does not necessarily involve Landau levels - for 
example, the in-plane magnetic held in graphene can induce excitonic instability 


(52). 


In a relatively weak magnetic held < eB < 4:irF% (F n is the pion decay 
constant), the correction to the chiral condensate can be evaluated using the 
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chiral perturbation theory (53) 


. f eB in 2 ( ( eB) 2 

S(-B) = S(0) 1+ 77 „ +0 1 


( 11 ) 


x (47 \ F% ))' 
this result demonstrates an increase of the condensate in magnetic field. Note 
that the dependence on the field in Eq. © is linear, and not quadratic as 
can be naively expected from power counting. This is because the pion loop 
diagrams describing the vacuum energy density in magnetic field yield logarithms 
involving the infrared cutoff that is usually provided by the pion mass. The chiral 
condensate can be obtained from the energy density by differentiating w.r.t. the 
quark mass, or, using the GMOR relation, w.r.t. the pion mass squared. As 
a result, the chiral condensate becomes inversely proportional to the infrared 
cutoff. However in an external magnetic field, the appropriate physical infrared 
cutoff becomes eB > - so the naive quadratic dependence on eB gets replaced 


by the linear one (53). In the absence of Goldstone modes, the dependence of 


the condensate on eB is quadratic (54), as naively expected. This discussion 
illustrates why an external magnetic field has a big effect on QCD interactions 
- it provides a physical (and anisotropic) infrared cutoff, and QCD possesses a 
large sensitivity to the infrared regularization. 


The lattice studies of QCD in magnetic field indeed indicate (55) the expected 


increase 


of the chiral condensate in magnetic field consistent with Eq. ©• On 
the basis of magnetic catalysis, it has also been widely expected that the restora¬ 
tion of chiral symmetry in QCD will be delayed to higher temperatures in an 
external magnetic field. However the first principle lattice QCD x QED results 
indicate the opposite trend, see Fig. [4]from (56). The lattice results show that, 
surprisingly, a phenomenon opposite to magnetic catalysis takes place in QCD 
matter - this has been called the “inverse magnetic catalysis”; the possible origins 
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of this effect have been discussed in (57,8). 

The properties of the QCD equation of state in a background magnetic field 


have recently been studied in (58), with intriguing results. In particular, while 


QCD matter appears weakly diamagnetic at low temperatures (consistent with 
the dominance of pions), around and above the transition temperature the re¬ 
sponse is paramagnetic (58). This means that the QCD plasma is susceptible to 
magnetization, e.g. in the form of helical magnetic configurations triggered by 


the chiral magnetic instability, see Section 7.1 The studies of QCD matter in 


magnetic field are discussed in much more detail in recent reviews (59 


3.2 Instantons in magnetic field 

Instantons (1) are configurations of the gluon field and thus do not interact 
directly with an external Abelian magnetic held. However quarks carry both 
electric and color charge and couple to both magnetic and instanton fields. A 
magnetic held introduces Landau levels in which the fermion zero modes have a 
dehnite spin projection on the magnetic held. In the instanton background the 
fermion spectrum also has zero modes, with chiralities determined locally by the 
local topological charge. 

The combination of background helds of the instanton and an Abelian mag¬ 
netic held thus leads to a surprisingly intricate and rich structure in the Dirac 


spectrum studied in (60). The result is driven by the competition of two effects: 


the spin projection produced by a magnetic held and the chirality projection pro¬ 
duced by an instanton held. Of particular interest is the emergence of electric 
dipole moment in a magnetic held that signals the violation of P and CP sym¬ 


metries localized around the instanton (60). The easiest way to understand this 
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phenomenon is to use the chiral representation of the quark spinors and to write 
down the electric dipole moment as 


°3 = -V’LC r 3'i/ , L + 


( 12 ) 


magnetic field is directed along X 3 , and a is the spin operator. As can be seen 
from Eq|12[ in the absence of chirality asymmetry the contributions from left- 
and right-handed quarks cancel each the electric dipole moment is equal to zero. 
However in the presence of the instanton the left-right symmetry is lost, and 
this leads to the emergence of the electric dipole moment. In the limit of strong 
magnetic field eB p(p is the instanton size), the Dirac spectrum can be 


evaluated exactly (60). 


Lattice calculations in QCD x QED confirm the emegence of electric dipole 
moment around the topological fluctuations of gluon field (6lJ 62,63 ,j64||65). Fig. 
[5] shows the electric dipole emerging around a topological fluctuation in QCD 


matter at temperature T = 113 MeV (65). 


3.3 Chern-Simons diffusion in magnetic field 


The Chern-Simons diffusion rate for strongly coupled J\T = 4 plasma in the pres¬ 


ence of an external magnetic field was computed via holography in Ref. (66) 


The dual gravity is characterized by a full solution of asymptotically AdS§, five¬ 
dimensional Einstein-Maxwell system with a constant magnetic flux. This solu¬ 
tion takes into account all of the back-reaction of the magnetic field on the metric 
and allows one to work with arbitrarily strong magnetic fields; it was introduced 


in (67), see also (68). 
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In a weak magnetic field, to leading order in B 2 one gets ( 66 ): 


r (9 2 N)\ 4 

T ° 8 = 


' 1 B 2 ,„(B^ 

1 + 6 T±T* + 0 {t* 


(13) 


note that at zero magnetic field this expression coincides with Eq]10| The pre¬ 
factor of B 2 /T 4 is numerically small ~ 2 ICC 3 , which suggests that the effect of 
magnetic fields on the sphaleron rate in heavy ion collisions may be neglected. 

In the strong magnetic field limit the diffusion rate is 


(n 2 IV Y 2 

T(B, T) = J - B T 2 


384 v/3- 


B » T 2 . 


7T 


(14) 


The result Eq[l4] can be understood as a consequence of dimensional reduction 
in strong magnetic field. In holographic approach, this is manifested through 


the emergence of a (2+l)-dimensional BTZ black hole in the bulk ( 66 ). It is 
interesting that this dimensional reduction applies even to the fields that are 
not charged under the external magnetic field - this is a consequence of strong 
coupling. 


4 The chiral magnetic effect (CME) 


The term Chiral Magnetic Effect (11,45) refers to the generation of electric cur¬ 
rent induced by the chirality imbalance in the presence of an external magnetic 
field: 


J= ^2 ^ B, (15) 

where ^5 is the chiral chemical potential. In QCD matter subjected to an external 
magnetic field or rotation, topological fluctuations thus induce fluctuations of the 
electric dipole moment (69 |70fllf 45l); we have already encountered a phenomenon 
of this type in section |3.2| addressing QCD instanton in a magnetic background. 
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The physics of CME and related phenomena has been recently reviewed in © 
which also contains an extensive set of references. In this paper we will thus limit 
ourselves only to a brief summary followed by the discussion of recent CME- 
related developments in heavy ion physics (section [5J, condensed matter physics 
(where we describe the first observation of CME in section [6]), physics of neutron 
stars (section [7.1[ ) and the Early Universe (section [7.2[ ). 

The CME is a macroscopic quantum effect - it is a manifestation of the chiral 
anomaly creating a collective motion in Dirac sea. Because the chirality imbalance 
is related to the global topology of gauge fields, the CME current is topologically 
protected (71) and thus non-dissipative even in the presence of strong interactions. 
As a result, the CME and related quantum phenomena affect the hydrodynamical 
and transport behavior of systems possessing chiral fermions, from the quark- 
gluon plasma to Dirac semimetals. 

The persistence of CME at strong coupling and small frequencies makes the 
hydrodynamical description of the effect possible, and indeed it arises naturally 


within the relativistic hydrodynamics as shown by Son and Surowka (72). The 
quantum anomalies in general have been found to modify hydrodynamics in a 


significant way, see (73) for a review. The principle of “no entropy production 
from anomalous terms” can be used to constrain the relativistic conformal hydro¬ 
dynamics at second order in the derivative expansion, where it allows to compute 


analytically 13 out of 18 anomalous transport coefficients (74). 

Anomalous hydrodynamics has been found to possess a novel gapless collective 


excitation - the “chiral magnetic wave” (75), see also (76). It is analogous to 
sound, but in strong magnetic field propagates along the direction of the field 


with the velocity of light (75). The chiral magnetic wave is the hydrodynamical 
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mechanism of transporting the CME current; it transforms an initial chiral or 
electric charge fluctuation into a macroscopic observable asymmetry in the dis¬ 


tribution of electric charge (77,78). Very recently, the chiral magnetic wave has 


been derived from the chiral kinetic theory (79). In particular, in the presence 


of chirality-flipping transitions, the chiral magnetic wave at frequencies smaller 
than the transition rate has been found to give rise to a diffusive vector mode 


(79). The chiral magnetic wave thus can provide an appropriate description for 


the phenomena discussed in section [6j 


5 CME in heavy ion collisions 

Relativistic heavy ion collisions produce hot QCD matter characterized by strong 
fluctuations of topological charge. In addition, the colliding ions generate strong 


magnetic fields eB ~ 0(10 m%) (11) that are directed, on the average, orthogo¬ 
nally to the reaction plane, see Fig. [6j 

Let us assume that the produced QCD matter contains a chiral charge produced 
e.g. due to a local imbalance between the sphaleron and anti-sphaleron transitions 
discussed in section [3| According to Eq. 15, there will then be a fluctuation of 
electric current that will result in a charge asymmetry relative to the reaction 
plane. This asymmetry will lead to the presence of P-odd harmonics in the 
angular distribution of the produced positive and negative hadrons: 

dN± 


dcj) 


1 ± 2asin(A(/>) + .... 


(16) 


with the asymmetry parameter a ~ 1% (69); Acf> is the azimuthal angle relative 


to the reaction plane, see Fig. [6] Since QCD does not violate P and CP globally, 
the averaged over many events value of the asymmetry is zero, (a) = 0. However 
the charge asymmetry fluctuates event-by-event, and the signature of CME is the 
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presence of dynamical fluctuations that exceed the statistical ones. 

The number of charged hadron tracks in a single event (although sufficient to 
determine the reaction plane) is not large enough to allow a statistically sound 
extraction of the coefficients a±, so one has to sum over many events. Because 
the sign of the charge asymmetry should fluctuate event by event, it was proposed 


by Voloshin (80) to measure the variable 


(cos (cj) a + 4>/3 — 2 'f’Rp)) = (cos A(f> a cos A(j)p) — (sin A 4> a sin A, (17) 

where the indices a and /3 denote the charge of hadrons, and T/j p is the reaction 
plane angle, see Fig. | 6 j The quantity Eq. 17 can be measured with a very high 
precision, and is directly sensitive to the presence of parity-odd fluctuations; it 
has an added benefit of not being sensitive to the reaction plane-independent 
backgrounds. There is a price to pay however - the observable itself is P-even, 
and so one has to carefully examine the possible backgrounds; see recent review 


(81) for a detailed discussion. 


In 2009, the STAR Collaboration at Relativistic Heavy Ion Collider at BNL 


reported (82,83) the measurement of the fluctuation of charge asymmetry Eq. 17 


at y/swN = 200 GeV; the order of magnitide of the measured asymmetry was con¬ 


sistent with expectations (69). Recently, the ALICE Collaboration Large Hadron 


Collider at CERN has reported (84) the measurement of the charge asymmetry 
at y/SNN = 2.76 TeV, see Fig. [7j The magnitude of the observed asymmetry is 
close to the one previously observed at RHIC - this is consistent with the fact 
that magnetic flux through the QCD plasma is approximately energy-independent 
(magnetic field is proportional to the 7 factor, and the longitudinal extension of 
the produced matter at early time is Vt)- 

The backgrounds to the measured charge cumulants from local hadron corre- 
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lations are expected to contribute to both the second harmonic Eq{T7] and the 
higher harmonics. Fig. [7j right shows the third harmonic as measured by the 
ALICE Collaboration - one can see that it is quite small, about an order of 
magnitude smaller than the second harmonic that is sensitive to CME. 

The propagation of the chiral magnetic wave in QCD fluid requires the restora¬ 
tion of chiral symmetry ( 75fl77 ) - therefore one expects a reduction in the strength 
of charge asymmetries when the collision energy decreases below the energy 
needed to induce the QCD transition. This is consistent with the results of 
STAR Collaboration from the beam energy scan at RHIC (85), see Fig. [8j 

A fully quantitative theoretical approach to the description of charge asymme¬ 
tries requires the use of relativistic hydrodynamics that includes the terms arising 
from the chiral anomaly, supplemented by the initial conditions describing topo¬ 
logical fluctuations at the early stage of a heavy collision. The study of that 


kind has recently been performed in (78); the initial condition is provided by the 


fluctuating longitudinal “glasma” fields with E a B a / 0 (86 87 88). The snapshot 


of the resulting chiral and electric charge densities in the QCD fluid is shown in 
Fig- § The computed (with and without anomalous terms) charge asymmetries 


are shown in Fig. 10 in comparison to the STAR experimental results (82). One 


can see that i) chiral anomaly has a big effect on charge asymmetries and ii) 
the results of the computation including the anomaly (and thus the CME) agree 
with the data within a factor of two, possibly leaving room for some background 
contributions. 

One of the remaining sources of uncertainty in theoretical calculations is the 


duration of magnetic held in QCD fluid, see (12) for a review. However this 


uncertainty can be reduced by the study of directed how of charged hadrons 
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away from mid-rapidity, as proposed recently in (89). The magnetic field is 


also expected to contribute to the photon and dilepton production through the 
“magneto-sono-luminescence”: the conversion of phonons into real or virtual pho¬ 


tons in a magnetic background (90,91) 


6 CME in condensed matter: the first observation 


The recent discovery (92 93 94) of Dirac semimetals with chiral quasi-particles 


opens a fascinating possibility to study the chiral magnetic effect (CME) in con¬ 
densed matter experiments. Very recently, the first observation of CME was 
reported through the measurement of magneto-transport in zirconium pentatel- 


luride, ZrTe .5 (14). 


The signature of the CME in Dirac systems in parallel electric and magnetic 
fields is the positive contribution to the conductivity that has a quadratic depen¬ 
dence on magnetic field (|45|95,96). This is because the CME current Eq. 15 


is 


proportional to the product of chiral chemical potential and the magnetic field, 
and the chirality imbalance in Dirac systems is generated dynamically through 
the chiral anomaly with a rate that is proportional to the product of electric and 
magnetic fields. As a result, the longitudinal magnetoresistance becomes negative 


(95,96), with a definite telltale dependence on magnetic field. In the framework 


of holography, the negative magnetoresistance was recently reproduced in (97). 


Let us explain how this mechanism works in Dirac semimetals in more detail. 

In the absence of external fields, each Dirac point contains left- and right- 
handed fermions with equal chemical potentials, hl = M R = 0. If the energy 


degeneracy between the left- and right-handed fermions gets broken, we can pa¬ 
rameterize it by introducing the chiral chemical potential ns = ([Ir — Hl)/ 2. The 
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corresponding density of chiral charge is then given by (45) 


P5 = 


4 


37T 2 t> 3 2>v 3 


+ T 2 + ^ 


7T 


(18) 


The chiral anomaly in parallel external electric and magnetic fields generates 
the chiral charge with the rate given by 

dp 5 < 


dt AiT 2 h 2 i 


-E-B. 


(19) 


The left- and right-handed fermions in Dirac semimetals can mix through chirality¬ 
changing scattering, and this process depletes the maximal amount of chiral 
charge that can be produced. Denoting the chirality-changing scattering time 
by Ty, we thus get the equation 


dp 5 _ _' 

dt 4ir 2 ti 2 c 


E-B- 


P5 

T V ' 


( 20 ) 


The solution of equation (20) at f ry is 


e 

47 T 2 h 1 c 


P5 = , 0 , 2 E ■ B T V . 


( 21 ) 


According to (18), this leads to a non-zero chiral chemical potential p§ = (pr — 
Pl)/2 (assuming that p§ <C p,T): 


3 v 3 e 2 E-B 


P5 ~ 


4vr 2 h 2 c T 2 + 4. TV ' 


( 22 ) 


The formulae (15) and (22) yield the expression for the CME current: 


ji _ 6 3 e V Ty j f, — jfc pfc 

J CME - —T oi“ Tf . „2 & t, = cr CME . 


irh 8 he 7T 3 X 2 _|_ p 


2 


(23) 


We see that the CME is described by the conductivity tensor <7y ME ~ B l B k . 
When the electric and magnetic fields are parallel, the CME conductivity is 

e 2 3 e 2 v 3 


°CME ~ 


T V 


7T h 8 he 7T 3 X 2 + —■ 


B . 


(24) 
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Since the CME current is directed along the electric field, it will affect the mea¬ 
sured conductivity - the total current will be the sum of the Ohmic and CME 
ones: 

J = <A)hm + JcME = (<7ohm + ^CMe) E, (25) 


where ff C ME = ^cme- 

If the electric and magnetic fields are parallel, there is no conventional contribu¬ 
tion to magnetoresistance induced by the Lorentz force. The magnetoresistance 


p = 1/a following from (25) is negative, with a characteristic quadratic depen¬ 
dence on magnetic field. It is this negative MR with a quadratic dependence on 


magnetic field that was observed in ZrTes (14). 

ZrTes is a layered material that crystallizes in the layered orthorhombic crystal 
structure. Prismatic ZrTe6 chains run along the crystallographic a-axis and are 
linked along the c-axis via zigzag chains of Te atoms forming two-dimensional 
(2D) layers stacked along the 6-axis into a crystal. 

The requirement for observation of the CME is that a material has a 3D Dirac 
semimetal-like (zero gap), or semiconductor-like (non-zero gap) electronic struc¬ 


ture. Figure 11, left shows angle-resolved photoemission spectroscopy (ARPES) 


data (14) from a freshly cleaved (o — c plane) ZrTes sample. The states form¬ 
ing the Fermi surface disperse linearly over a large energy range, both along the 
chain direction (Figure [Tl] left) and perpendicular to it (14), indicating a Dirac- 
like dynamics of carriers for the in-plane propagation. The velocity (the slope 
of dispersion) is very large in both the chain direction, v a — 6.4 eVA(~ c/300), 
and perpendicular to it, v c — 4.5 eVA- these values are close to the velocity in 
graphene, v — c/300. 


However, the in-plane electronic structure of ZrTes cannot be described by 
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a single (anisotropic) cone, especially away from very low energies. The states 
are doubled - this is because the crystal contains two layers per unit cell - and 
the simplest description is that the bi-layer splitting creates two cones (bonding¬ 
antibonding), separated in energy by ~ 300 meV possibly with a small gap at 


the degeneracies, as shown in Figure llll right, from (14). 


A negative magnetoresistance (MR) is observed (14) in ZrTes for T < 100 K, 
increasing in magnitude as temperature decreases. The magnetic field dependence 
of the negative MR can be nicely described by the CME contribution to the 


electrical conductivity, given by (24) - this is illustrated in Fig. 12, right for 


T = 20 K, from (14). At 4 Tesla, the CME conductivity is about the same as 


the zero-field conductivity. At 9T, the CME contribution increases by ~ 400%, 
resulting in a negative MR that is much stronger than any conventional one 
reported at an equivalent magnetic field in a non-magnetic material. Inset in 


Fig. 12 shows the temperature dependence of the pre-factor of B 2 which decreases 
with temperature faster than 1/T, consistent with the CME. 

At very low field, the data show a small cusp-like feature. The origin of this 
feature is not completely understood, but it probably indicates some form of anti¬ 
localization. Similar feature has also been seen in the Dirac semimetal CdsAs 2 


(98), where an indication of negative magnetoresistance at larger fields has also 


been reported. It would be interesting to understand the origin of this positive 
MR at weak fields. 


As pointed out in (14), the studies of CME can be extended to a broad range of 


materials since 3D Dirac semimetals often emerge at quantum transitions between 
normal and topological insulators, including topological crystalline insulators. 
Dirac semimetals also open new possibilities for photonics since they are expected 
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to possess well pronounced plasmon excitations in the THz frequency range, with 


universal properties (99) 


7 Topology and the Universe 


7.1 CME in neutron stars 

Supernova collapse is known to provide a kick to the emerging proto-neutron 
stars. The velocity of the produced neutron stars reaches v = 1083^go 3 km/s, 
as measured by Very Long Baseline Array (100) for the pulsar B1508+55. The 
mechanism of this kick remains a mystery. Several of the proposed explana¬ 
tions involve the interplay of neutrino physics with the magnetic field, e.g. the 


anisotropy of neutrino emission (101) or neutrino oscillations biased by magnetic 


fields (102). 


It has been proposed in (103) that the neutron star kick can be provided by 


the chiral magnetic current arising from the asymmetry between left- and right- 
handed electrons due to the parity violation in weak interactions. The authors 
estimate that CME throughout the typical age t ~ 40,000 years of a neutron star 
can provide it a velocity on the order of v ~ 1000 km/s. 

A different CME-related mechanism of the neutron star kick was proposed in 
(104). The idea of (104) is based on the “chiral separation effect” that is dual 
to the CME - it is a flow of the axial current induced by the baryon chemical 
potential. Because neutrino scattering can tell the difference between left- and 
right-handed fermions, this flow of axial current induces an anisotropy in neutrino 
emission and thus the desired kick of the neutron star. 

It would be interesting to explore the propagation of the axial current using 
the anomalous hydrodynamics, in analogy to previous studies of the “chiral mag- 
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netic wave” (75.77) in heavy ion collisions. The first steps in using anomalous 


hydrodynamics have already been made (105). In addition to magnetic field, the 


rotation of neutron stars also has to be taken into account since they bring in the 


effects related to gravitational anomaly (106) 


Another long-standing problem in the physics of neutron stars is the origin of 
their strong magnetic fields. The Maxwell-Chern-Simons theory that provides 
(71) an effective low energy theory of CME in the presence of chiral chemical 


potential ^5 displays a “chiral magnetic instability” (107.108,109 110) of the 
magnetic field at small momenta k < e 2 This instability leads to the 
possibility of producing helical magnetic fields characterized by non-zero value of 
magnetic helicity (Chern-Simons 3-form) f d 3 x A- B. During the collapse of neu¬ 
tron stars, left-handed electrons are captured by protons due to weak interactions. 
This induces a chiral imbalance that through the chiral magnetic instability can 


power the growth of a helical magnetic field, as proposed in ( 111 ) 


A more detailed investigation (112) however revealed the important role of 


electron mass that dampens the chiral magnetic instability, at least in the for¬ 


mulation of (111). It remains to be seen if the chiral magnetic instability can be 


saved in a more elaborate treatment taking into account the energy balance in 
magnetic and gravitational fields of the staiQ 


7.2 CME and the primordial magnetic fields 


An intriguing application of CME that stimulated some of the early work (113 


109,114,110) as described in review ©> is the generation of primordial magnetic 
fields in the Early Universe. The chirality imbalance in the primordial electroweak 


1 I thank S. Reddy for a stimulating discussion of this issue. 
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plasma can be converted by the chiral anomaly into a helical magnetic field 
configuration with non-zero Chern-Simons number. 

To explain the presence of relatively strong magnetic fields and magnetic he- 
licity in the present day Universe, one has to understand the transfer of magnetic 


helicity from small to large scales (115), study of the conversion of the electroweak 
plasma into a horizon-scale helical magnetic held (116), and the realization that 
leptogenesis can give rise to right-handed helical magnetic held that is coherent 


on astrophysical length scales (117). Similar ideas have been developed in the 


series of papers (118 119,120). The chiral vortical effect has been found to lead to 
the production of the helical magnetic held in the turbulent electroweak plasma 
(121), whereas the CME amplihes the growth of the held. 

A particularly intriguing option for the Universe is the inhationary expansion 
driven by a pseudo-scalar inhaton - for example, a pseudo-Nambu-Goldstone 


boson (122) or an axion (123). In this case the entire Universe would be in a 
parity-odd state, and the coupling to gauge helds would induce CME currents 
on cosmological scales. In particular, this cosmological parity violation would be 
imprinted in the Cosmic Microwave Background, as discussed recently in (124j 


125 

126 

127 

128 

129 

130) 


8 Outlook 

Topological held conhgurations dehne many fundamental properties of the phys¬ 
ical world, including spontaneous breaking of chiral symmetry and confinement 
in QCD. As argued throughout this review, coherent magnetic helds provide an 
ideal probe of QCD topology. For the investigation of non-perturbative dynamics 
of QCD, external magnetic helds may well play a role similar to the one played by 
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deep-inelastic scattering in establishing the short distance perturbative behavior 
of quarks and gluons. 

Strong magnetic fields are accessible in heavy ion collisions and in numerical 
simulations on the lattice; this opens unique possibilities for testing theoretical 
ideas about QCD topology and gaining a deeper understanding of strong inter¬ 
actions. The presence of strong magnetic fields (up to ~ 10 15 G) in neutron stars 
provides an access to magnetized cold dense QCD matter. The knowledge gained 
about topological effects in non-Abelian plasmas can advance the understanding 
of the baryon asymmetry and magnetic helicity in the Universe. 

Finally, the recent discovery of Dirac semimetals and the observation of the 
chiral magnetic effect in ZrTes open an exciting possibility to study 3D topological 
field-theoretic phenomena in condensed matter systems, with a potential for a 
wide array of applications. 

I thank M. Chernodub, M. D’Elia, G. Endrodi, Z. Xu and the STAR Collabo¬ 
ration for the permission to use graphics from their papers and helpful comments 
on the manuscript. This work was supported in part by the U.S. Department of 
Energy under Contracts DE-FG-88ER40388 and de-sc0012704. 
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Figure 1: The temperature dependence t = T/T c — 1 of the ratio x(T)/%(0) of 
topological susceptibilities at non-zero and zero temperature T in SU(N) gauge 
theory; from @ 0 - 
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Figure 2: The temperature dependence t = T/T c — 1 of the 62 moment of the 
topological charge distribution in SU(N) gauge theory for N = 3 and N = 6 ; from 


( 16 ). 
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Figure 3: The conjectured phase diagram of QCD in the temperature - chiral 


chemical potential plane; from (47). 
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Figure 4: The phase diagram of QCD in the temperature - magnetic field (T, eB) 


plane determined from Nf = 2 + 1 lattice QCD; from (56). 
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Figure 5: The electric dipole structure associated with the topological charge in 
QCD matter (left panel); the spatial correlation of the electric and topological 
charges along the direction of magnetic held at temperature T = 113 MeV; from 


( 65 ). 
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the reaction plane and of a produced hadron; from (85). 
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Figure 7: Charged hadron correlator as a function of collision centrality measured 


by the ALICE Collaboration; from (84). Left panel: the second harmonic Eq. 


17 right panel: the third harmonic of the azimuthal angle distribution. 
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Figure 8: The three-point correlator measured by STAR Collaboration in Au-Au 
collisions at RHIC as a function of centrality at different collision energies; from 


( 85 ). 
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Figure 9: Distributions of the chiral (left) and electric (right) charge densities 
in the transverse plane at mid-rapidity and a proper time r = 1.5fm of a AuAu 
heavy ion collision at yfs = 200 GeV per nucleon as computed in anomalous 


hydrodynamics; from (78). 
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Figure 10: The three-point correlator 7 for same and opposite charge hadrons as 


computed using anomalous hydrodynamics compared to the results from STAR 


Collaboration in Au-Au collisions at y/s = 200 GeV per nucleon; from (78). 
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Figure 11: Electronic structure of ZrTes: the valence band dispersion along the 
chain direction measured by Angle Resolved Photo Emission Spectroscopy (left); 


schematic view of the in-plane low energy electronic structure; from (14). 
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Figure 12: Longitudinal magnetoresistance of ZrTes: the experimental measure¬ 
ment (open red points) compared to the prediction from the chiral magnetic effect 


(blue line); from (14). 






